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An approximate analytical scheme of the dynamical mean field theory (DMFT) is developed for 
the description of the electron (ion) lattice systems with Hubbard correlations within the asymmet- 
ric Hubbard model where the chemical potentials and electron transfer parameters depend on an 
electron spin (a sort of ions). Considering a complexity of the problem we test the approximation in 
the limiting case of the infinite-?/ spinless Falicov-Kimball model. Despite the fact that the Falicov- 
Kimball model can be solved exactly within DMFT, the densities of states of localized particles have 
not been completely investigated off half filling. We use the approximation to obtain the spectra of 
localized particles for various particle concentrations (chemical potentials) and temperatures. The 
effect of a phase separation phenomenon on the spectral function is considered. 
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I. INTRODUCTION 

The Falicov-Kimball modeli was introduced to de- 
scribe the thermodynamics of metal-insulator transitions 
in compounds that contained both itinerant and localized 
quasiparticles. The spinless Falicov-Kimball model is the 
simplest example of an interacting fermionic system that 
displays numerous phase transitions. Despite its rela- 
tive simplicity, the analysis of the model is very complex 
and still a lot of open problems remain. The model can 
describe many phenomena such as metal-insulator tran- 
sition, ferromagnetism, antiferromagnetism, phase sepa- 
rations, etc. 

In the dynamical mean-field theory 2,3 which is exact 
in infinite dimensions^ the Falicov-Kimball model can 
be solved exactly by calculating the grand canonical po- 
tential and the single-particle Green's function^Li. The 
thermodynamics of the Falicov-Kimball model is well in- 
vestigated. The phase transitions between homogeneous 
phases or phase separations are investigated using the 
mentioned exact methods^ as well as using the effec- 
tive Hamiltonian at the large-C/ limi1*i& (see Ref. ^] for 
detailed review of the Falicov-Kimball model). 

The problem of an evaluation of the spectral function 
of localized particles is more complex than the investiga- 
tion of thermodynamics and the phase transitions. The 
/ spectral function related to the localized particles was 
calculated at half-filling more than ten years agoiS. Re- 
cently, the exact scheme has been extended to the case 
of different particle concentrations^ 3 -, but due to com- 
putational difficulties the densities of states of localized 
particles have not been completely investigated off half 
filling. So, we suggest here an approximate analytical 
scheme within the dynamical mean-field theory for cal- 
culating the / spectral function at different particle con- 
centrations and temperatures. 

We consider the asymmetric Hubbard model describ- 
ing the dynamics of two types of particles (ions, elec- 
trons or quasiparticles) as a generalization of the Falicov- 
Kimball model. The Hamiltonian of the asymmetric 



Hubbard model in a second quantization has the follow- 
ing form 

H = - ^n la + U n *T n a + % a \<T a ia > W 

ia i ija 

where rii a — a\ a a ia and the motion of particles is de- 
scribed by the creation (a icr ) and annihilation (a irT ) op- 
erators. The chemical potentials \i a and the transfer pa- 
rameters tfj depend on a sort of particles (an electron 
spin). The value U describes the local on-site repulsion. 

The asymmetric Hubbard model was proposed for the 
description of mixed- valence compounds 1 —. This model 
can also be used for the investigation of the lattice sys- 
tems having an ionic conductivity with two types of ions. 
When a single lattice site can be occupied only by one 
ion, the nontrivial limit U — > oo should be considered. 
Is this case, various thermodynamic regimes can be real- 
ized. The chemical potentials or concentrations of parti- 
cles of different sorts can be fixed independently. In this 
context the model can be investigated in the presence of 
the external field corresponding to the difference between 
the chemical potentials of different sorts. 

A number of methods for describing the strongly cor- 
related electron systems has been developed within the 
dynamical mean field theory. However, all these meth- 
ods have various restrictions. The quantum Monte-Carlo 
methodi£*i£ii£ is numerically exact but has severe prob- 
lems at low temperatures and for high repulsion strength 
U. The exact diagonalization mcthodi2*i£ is restricted to 
a small number of orbitals. Among the numerical tech- 
niques the most reliable one at low temperatures is the 
numerical renormalization group method 20 . For example 
its extension was used for the description of the ground 
state of the standard Hubbard modelSi. 

Besides the numerical approaches the development 
of analytical approximations for the infinite-dimensional 
model still remains necessary. The Hubbard model in 
the large-C/ limit was investigated using the non-crossing 
approximation^ 2 ^. Many approximations were devel- 
oped for the weak-coupling regime, for example, the 
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Edwards-Hertz approach 24 i 25 i 26 i 27 . The alloy-analogy 
based approximations 2 ^ 2 ^*^ do not take into account the 
effect of scattering processes on forming the energy band 
and cannot be used for the investigation of spectra of 
the asymmetric Hubbard model. In the Falicov-Kimball 
limit, the alloy-analogy (A A) and modified alloy- analogy 
(MAA) approximations give the density of states of lo- 
calized particles in the form of a delta function. The 
scattering processes should be taken into account for a 
correct description of the broadening of this peak. For 
example, it was the Hubbard-III approximation^ that 
originally included the electron scattering for the half- 
filled Hubbard model into the theory. 

We use and improve the approximate analytical ap- 
proach originally proposed for the Hubbard models and 
extended to the asymmetric Hubbard model^. In this 
method the single-site problem is formulated in terms 
of the auxiliary Fermi-field. The approach is based on 
the equations of motion and on the irreducible Green's 
function technique with projecting on the basis of Fermi- 
operators. This approach gives DMFT equations in the 
approximation which is a generalization of Hubbard-III 
approximation and includes as simple specific cases the 
AA and MAA approximations. 

The approximation is tested on the infinite-C/ spinlcss 
Falicov-Kimball model. We use the approximation to ob- 
tain the densities of states of localized particles for var- 
ious chemical potentials (concentrations) and tempera- 
tures. The dependences of the chemical potentials on 
the particle concentrations are calculated using the densi- 
ties of states as well as thermodynamically by calculating 
the grand canonical potential^. To calculate the particle 
spectra at low temperatures the phase separation should 
be taken into account. 

In Section II, we review the formalism of DMFT with 
the use of the auxiliary Fermi-field and the approxi- 
mate analytical scheme based on the projecting tech- 
nique and the different-time decoupling procedure. In 
Section III and in Appendix, the exact relations between 
some Green's functions are derived, which allows us to 
find the projecting coefficients using only the single par- 
ticle Green's function and the coherent potential. The 
results are discussed in Section IV, followed by our con- 
clusions in Section V. 



with the coherent potential J a (r — r') which has to be 
self-consistently determined from the conditions 



G a (uj„) 
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G -K> k ) = g -ir„\ 

G a K) = G^K) = ^ E G<T (^«> k )> (5) 

k 

where G a is the one-particle Green's function, S CT is the 
total irreducible part which does not depend on the wave 
vector k. The sum over k in Eq. is calculated by the 
integration with the density of states (DOS) (a Gaussian 
DOS for an infinite-dimensional hypercubic lattice and a 
semielliptic DOS for a d = oo Bethe lattice). 

It was shown in Refill that the single-site problem can 
be formulated in terms of the auxiliary Fermi-operators 
£ CT describing the creation and the annihilation of par- 
ticles in the effective environment. The problem is de- 
scribed by the following Hamiltonian 



H cS = H Q + J2V 



(6) 



The approach does not require an explicit form of the 
environment Hamiltonian H^. The environment is de- 
scribed by the coherent potential given as the Green's 
function for the auxiliary Fermi-field with the unper- 
turbed Hamiltonian H^: 



J CT (u,)=27«<Cia» 



t\\€ 



(7) 



The particle creation and annihilation operators are 
expressed in terms of Hubbard operators 



a a = X w + (X c 

on the basis of single-site states \nA,ns) 

|0) = |0,0), \A) = \1,0), 

|2) = |1,1), |B) = |0,1>, 



(8) 



(9) 



where the following notations for sort indices are used: 
er = B, £ = + for cr = A and a = A, ( — — for a = B. 
In this representation the local Hamiltonian Hq of the 
asymmetric Hubbard model is 



H 



[/i ff (X CTCT + X 22 )] + UX 



22 



(10) 



II. FORMALISM 



In the dynamical mean-field theory the infinite-dimen- 
sional lattice model is mapped on the single-site problem 



e -0H _^ e -f3H clf = e -0H o 



x Texp 



dr /VVj ff (r - r')al{r)a a {r') 



(2) 



and the two-time Green's function G a (u>) = 2-K{{a a \a\,)) u 
is written as: 

G a = 2n[((X°°\X°°)) UJ +(((X °\X 2 *)) UJ 

+ Q{{X S2 \X M )) U + {{X^X 25 ))^] . (11) 

The Green's functions in Eq. ^2 are calculated using 
the equations of motion for Hubbard operators: 



(12) 
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The commutators (11211 are projected on the subspace 
formed by operators X 0(T and X a2 : 



eft 



aJX 0a + a 7 2 X a2 + Z~< . (13) 



The operators Z 0a ( a2 ' > are defined as orthogonal to the 
operators from the basic subspaceS^S: 



({Z 



0a(a2) j£-0o-(o-2) 



) }) = o. 



(14) 



These equations determine the projecting coefficients 
a 0cr(a2) w j 1 j c j 1 are expressed in terms of the mean value 



(15) 



Using this procedure by differentiating both with re- 
spect to the left and to the right time arguments, we come 
to the relations between the components of the Green's 
function G a and scattering matrix P a . In a matrix rep- 
resentation, we have 



G(j — Gq -\- Gq P(jGq , 



where 



f -o n f {{X°°\X 2 *)) 
Ua ~ \ {(X 92 \X^)) ((X 92 \X 2 °)) 

and nonperturbed Green's function GZ is 



(16) 



(17) 



D„ 



o a 29 r 



where 



(18) 



A pq = (X** + Xii), A Qa = l- n ff , A 2S = ng, (19) 



D„ 



{lo - a a ){uj - b a ) 



V? 



Ao a A 2 g 



(20) 



Clcr = ~fJ>a + -T-tPvi b a =U- flcr+ ^-fa- (21) 

The scattering matrix 



Prr = 27T 







A' 1 

Oct 

A~ 



{{Z°°\Z° )) {{Z 0a \Z 2s )) 
{{Z S2 \Z° )) {{Z S2 \Z 2S )) 



A 



Oct 





A' 1 



(22) 



being expressed in terms of irreducible Green's functions 
contains the scattering corrections of the second and the 
higher orders in powers of V a . The separation of the 
irreducible parts in P a enables us to obtain the mass 
operator M a = -P CT |ir and the single-site Green's function 
expressed as a solution of the Dyson equation 



G„ 



(1 - G^M^Gl 



(23) 



We restrict ourselves to the simple approximation in 
calculating the mass operator P a , taking into account 
the scattering processes of the second order in V a . In 
this case 



M„=P<?\ 



(24) 



where the irreducible Green's functions are calculated 
without allowance for correlation between electron transi- 
tion on the given site and environment. It corresponds to 
the procedure of different-time decoupling, which means 
in our case an independent averaging of the products of 
X and £ operators. 

Let us illustrate this approximation on the example of 
calculation of the following irreducible Green's function: 

/M = «(x°° + xnatu* 00 + x aa m • (25) 

According to the spectral theorem, this Green's function 
is related to the corresponding time correlation function, 
and according to the different-time decoupling we have 

(£(t)(x 00 + x™) t (x 0Q + x™)U" 

« ((X™ + X™ ) t (X 00 +X™))(£(t)U- (26) 

Calculation of these correlation functions in a zero ap- 
proximation 



{{X 00 + X aa ) t {X 00 + X aa )) 
» ((X 00 + X aa ) 2 ) = A 



Oct 



leads to the result 



i(u)=A 0a (ua))i = 



_Aqo_ 

2ttV : 



(27) 



(28) 



Using the above procedure we can obtain the final ex- 
pressions for the mass operator and the total irreducible 
part: 



A 0a 



A 2 g 



<■*> + Mcr — &<r{uj) <jJ + — U — ^ct(w) 



(29) 



where 



and 



£l„(u;) = J a [u) 



R a {u) = - 



(X° 



X 6 



Ao a A 2 g Ao a A 2 g 



-J S (UJ + fl a - fa) 



(30) 



{X™ - x ss 

2^ 

(X 00 + X 22 ) 



(X m - X 
2tt 



22 \ 



2 

du/Im J s (w' + ie) , QJ 
tanh 

UJ — Ld' — fig + fl a 2 

— 00 

Jg(U - fl a - Ha - t0) 
+ 00 

da/ImJ g (-q/-i £ ) tanh ^/ 

LU - LO' + fig + fl a — U 2 
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The approach used here for the approximate solution 
of the single-site problem can be called as the general- 
ized Hubbard-III (GH3) approximation. It becomes the 
standard Hubbard-III approximation in the case of the 
usual Hubbard model at half-filling with spin degenera- 
tion (n T =n L = 1/2, MT = p L = U/2, {X 00 ) = (X 22 ), 

{x aa ) = {x ss )y. 



s- 1 H= 



1/2 



1/2 



to 



U/2- 



3J a {uj) LO-U/2-3J a (io) 



(32) 



The function R a {u) describes band forming for parti- 
cles of sort a by the motion of particles of another sort 
a (scattering processes). The neglect of this contribu- 
tion (R a (oj) = 0) gives MAA approximation. If we put 
R a (uj) = and (p a = 0, the system is described within 
the simple AA approximation. 

In the limit of infinite repulsion U the following solu- 
tion of the single-site problem is obtained 8 : 



1 — n d 



(33) 



J S (iU + p a - (Igr) 



2?r 



ImJ,(q/ + i £ )da/ tanh ^ 

LU — LU 1 - p s + Pa 2 



The constant ip a = (£ 5 A ff0 ) + C,{X a2 S} s ) can be calcu- 
lated using the exact relation given in the next section. 
The average particle concentrations are calculated using 
the imaginary part of the Green's functions (density of 
states): 



p a (oj) = lmG a {uj + ie) 



(35) 



The self-consistency conditions (a set of equations © 
(0) relate the coherent potential J a to the Green's 

function G a . For the Bethe lattice with a semielliptic 

density of states 



we have 



W 2 

J a (u,) = ^G„(cu). 



(36) 



(37) 



In the case of the Falicov-Kimball model the unper- 
turbed bandwidth is zero (2Wb = 0, Jb(cj) = 0) for lo- 
calized particles, and the approach gives the exact equa- 
tion for the Green's function of itinerant particles Ga{w)- 
The density of state Pa(u) on the Bethe lattice is nonzero 
for \tu + Ha\ < WaVI n B- 



PaM = ^2 ^Wl{l-n B )-{u + p A y. (38) 



In this case the equations l|33(l . I|34|l give an explicit ap- 
proximate expression for the Green's function of localized 
particles: 

1 - riA 



ra 1-riA 1-riA 



(39) 



where 



d ( \ n B +n A , . 

Rb{u) = -z Ja(U + PB~ PA) 



Ub-Ua f ImJ A (^ + i £ )d^ tanh ^/^ 
2ir J lu - u)' - p A + pb 2 



III. THE PROJECTING COEFFICIENTS IN 
EQUATIONS OF MOTION FOR HUBBARD 
OPERATORS 

The average values (X pq ^ IJ ) are calculated using cor- 
responding Green's functions according to the spectral 
theorem 



dw 



s /3w + i 



2lm((UXn) 



UJ + i£ 



(41) 



These Green's functions can be calculated using the exact 
relation (|A.3(I derived in Appendix 



(42) 



For the asymmetric Hubbard model we have to calcu- 
late the coefficients V a <p a : 

Lp a = + C(* ff2 ^> = (6(A ff0 + cx 2a )) (43) 

which are expressed as 

'-+ 00 dLU 



2 Puj + I 



(44) 



p*{uo) = 2J a {u)((X° s - CX- 2 |A M + CA 2 ^» W 

= 2J ff ( W ) [((A 0ff \X*°)) U - ((X^\X^)) U } .(45) 

In the limit of infinite on-site repulsion U the state with 
double occupation is excluded and we have 



p%{u) = -J a (oj)G a (uj) . 

7T 



(46) 



For the Falicov-Kimball model it is possible to calculate 
the exact Green's functions for itinerant particles and it 
gives the exact expression for Vb<Pb- 

Let us note that previously this parameter was calcu- 
lated approximately using the Green's functions obtained 
by means of the linearized equations of motion and ne- 
glecting the irreducible parts&22i2ii. In the model with 
infinite U this approximate value is 



-J a {lu) — 
it lu - 



n B 



PB 



(47) 
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1,0 
0,5 



W A =1 

1) T=0.02 

2) T=0.3 

3) T=0.5 




FIG. 1: The exact values of ip (solid lines) are compared 
with the approximate ones calculated using the relation I47H 
(dashed lines) for various temperatures and tib- One can 
see that the approximation is applicable only for very high 
temperatures. 



and it corresponds to the approximation of the series 
(|A.9|1 where only the first term (the zero approximation 
for Ga) is taken into account. 

There is shown a comparison between the approximate 
and the exact values of ips in Fig. ^ The improvement 
given by the exact relation (|46|l allows us to investigate 
the system at low temperatures where summing up the 
whole series (IA.9I) is essential. 



IV. RESULTS 

Dependences of chemical potentials ha and /is on the 
particle concentrations are calculated using correspond- 
ing densities of states. DOS of localized particles is ob- 
tained as an imaginary part of the Green's function 



1 



ImGfi (oj + ie) 



(48) 



This Green's function has the correct analytic properties 
(relations between imaginary and real parts) within the 
considered approximations (GH3, MAA, AA). The ap- 
proximate DOS always has the correct sign and the sum 
rule is fulfilled: the integral of DOS over all frequency 
is equal to unity for finite U and is equal to 1 — n s for 
U — > oo when the upper band tends to infinity. 

For the Falicov-Kimball model the density of states 
of localized particles can be calculated exactl}iii*i2ii£ ; 
but numerical results were obtained mostly at half-filling. 
The constant ips is zero at half-filling (ua = tib = 1/2) 
because of the particle-hole symmetry. In this case we 
have the simple approximate solution of the single-site 
problem 



G B {u) 



oj — 2Ja(oj) 



lo 2 - U 2 /4-2ujJa(uj)' 



(49) 



This result at half-filling is independent of temperature, 
but for high temperatures and large values of U the ap- 
proximate scheme reproduces the exact results (Fig. [2J - 




FIG. 2: DOS of localized particles for the Falicov-Kimball 
model on the hypercubic lattice with noninteracting DOS 



pT(e) = tt 



Solid line 



1/2 exp(-e 2 ) at half-filling (n A = n B = 1/2). 
our approximation; dashed line - exact result—. 



For simplicity we restrict our investigation to the 
Falicov-Kimball model with the infinite repulsion U on 
the Bethe lattice. In this case the model describes the 
system with an average site occupation no more than 
unity. There is a homogeneous state for temperatures 
larger than critical T c ss 0.060WOL. For lower tempera- 
tures (T < T c ) there are various types of phase transi- 
tions depending on a thermodynamic regime&Siifi. 

First, we consider a possibility of describing the phase 
transitions using the approximate equations. For this 
reason, the behavior of dependences of the chemical 
potential [iB on the concentration of localized parti- 
cles n B is investigated. The phase transition is in- 
dicated by the thermodynamic unstable region where 
d/iB/driB < 0. When the critical temperature is ap- 
proached from below, the dependence ^t_e("-s) becomes 
monotonic (d^B/dns > 0). 

In Fig.[3]the approximate curves /is^s) are compared 
with the exact results obtained thermodynamically 8,33 . 
The alloy-analogy based approximations give the density 
of states of localized particles in a form of a noninteract- 
ing delta-function, which is correct in the atomic limit 
Ia = 0. Thus, the MAA approximation can give reason- 
able results only when the chemical potential fiA < (i.e. 
the concentration of itinerant particles tends to zero). 
This approximation shows the presence of a phase transi- 
tion. However, it overestimates the critical temperature. 
In Fig. |3| the exact curve at T — 0.060 is already mono- 
tonic but the MAA approximation shows the region with 
unstable concentration values, i.e. there is still a phase 
transition. 

The GH3 approximation, unlike the MAA approxima- 
tion, incorporates the scattering processes forming the 
energy band of localized particles. This is crucial for the 
calculation of thermodynamic quantities. The approxi- 
mate curves coincide with the exact ones in a wide range 
of temperatures for the concentrations larger than some 
value which depends on ha- In Fig-EK^A = 0) one can 
see a good agreement with the exact result at fiA = for 
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FIG. 4: DOS of localized particles within the GH3 approximation for various fiA and temperatures. Wa ~ 1; Wb = 0; U = oo; 
1) T = 0.06, 2) T = 0.2. 



concentrations > 0.5. At temperatures lower than the 
critical one the approximation clearly indicates the phase 
transition and gives the correct value for T c )c 0.060. 

Spectra of localized particles are plotted in Fig. In 
this case temperatures are higher than the critical one, 
so a homogeneous state is stable. The parameter values 
are chosen so that the approximation gives the correct 
thermodynamic (Ib{tib) relations. The energy band of 
itinerant particles (A) depends only on a concentration 
riB and its width is 2Wa\/T- -~nB- However, the band 
of localized particles (B) is generated by scattering pro- 
cesses and its spectral shape depends on the concentra- 
tion riB , the chemical potential ha (or the corresponding 
concentration ha) and temperature. 

Depending on the values of the chemical potential of 
itinerant particles [iA (or ua) there are two limiting cases 
with different properties of the localized particle spec- 
trum. For very small ua (negative ha) the system is 
close by its behavior to the atomic limit, i.e., the spec- 
trum pb{u) is in the form of a delta peak. The sharp 
peak is slightly broadened by the scattering of itinerant 
particles (Fig. ^V, HA = —0.4). The contrary case with 
the nearly filled bands is when the total concentration of 
particles tends to unity, ha+ub — ► 1 (Fig. HP, ^A = 0.4). 
In this case the peak vanishes because the contribution of 



the coherent potential of mobile particles Ja(w) and the 
corresponding function Rb{u) (|40(l becomes larger and 
the simple pole in (|39J) disappears. The spectrum ps(w) 
corresponds by the form to the lower Hubbard subband 
(or the lower subband for the half-filled Falicov-Kimball 
model obtained in Ref. IT^) with the chemical potential 
located in a gap in the strong coupling limit. The inter- 
mediate case with the broad band superimposed by the 
sharp peak is shown in Fig. 03. 

States with concentration values in some region (where 
d/itr/dria < 0) are thermodynamically unstable at low 
temperatures. So, the presence of phase transitions 
should be taken into consideration. There are phase tran- 
sitions between homogeneous phases with a concentra- 
tion jump in the thermodynamic regime with the fixed 
chemical potentials (pa — const, fiB — const). In this 
case the density of states changes instantly as the concen- 
tration jumps. If one of the concentrations (ha or ng) is 
fixed the homogeneous state is unstable and the phase 
separation takes place. The phase transitions for the 
Falicov-Kimball model were investigated in many works 
(see Ref. ITlj) ; the regimes mentioned above were inves- 
tigated for the Falicov-Kimball model using the exact 
thermodynamic equations in Ref. 0- 

Let us consider the thermodynamic regime with fixed 
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FIG. 5: DOS of localized particles within the GH3 approx- 
imation (Wa = 1, Wb = 0). A) At temperatures below T c 
(T = 0.059) the state is separated into two phases with differ- 
ent particle concentrations (1,2 - the spectra of each compo- 
nent; 3 - the superposition of 1 and 2). B) DOS for various 
temperatures (1 - T = 0.059, 2 - T = 0.061, 3 - T = 0.2, 4 - 
T= 1.0). 



values of \ia and ns . In the homogeneous state the spec- 
tral function of itinerant particles pa depends on ng and 
[iA (Eq.[3SJ) and is independent of temperature. However, 
in the phase separated state the spectrum of the whole 
system can be considered as a superposition of spectra 
of each component. The homogeneous state is unstable 
for the concentration values jib £ (nBi,fiB2) and the 
system is separated into two different phases with the 
concentrations of B particles nsi ~ 0.52 and ns2 ~ 0.78 
at T — 0.059 Wa. The concentrations of the components 
(riBi, fiB-i) depend on temperature^. Thus, the spectra 
Pa and ps are temperature dependent in the segregated 
phase. In Fig. [5{A) DOS of localized particles is plot- 
ted at temperature T = 0.059Wa < T c . In Fig. \RB) 
the spectrum at various temperatures is compared. The 
bandwidth is larger in the segregated state than in the 
homogeneous state (T > T c ) and is temperature depen- 
dent. 



V. CONCLUSIONS 

The approximate analytic approach within DFMT for 
calculating the single-particle Green's functions of the 
asymmetric Hubbard model is developed and improved. 
This approximation allows one to investigate the model 
for various concentration values. The method is tested 
on the Falicov-Kimball model with the infinite on-site 
repulsion. It is shown that for high enough tempera- 
tures or large concentrations of localized particles the 
approximate approach reproduces exact values of chemi- 
cal potentials. The approximation can correctly indicate 
the instability of a homogenous state and the presence of 
phase transitions. 

The generalized Hubbard-III approximation (GH3) 
partially includes the scattering of particles into the the- 
ory and describes the formation of the band of localized 
particles. In the infinite-l/ limit the spectrum of localized 
particles is obtained for various particle concentrations 



and temperatures. The form of this spectrum continu- 
ously changes from a delta peak to the characteristic form 
of the lower subband of the spectrum in the Hubbad-III 
approximation when the chemical potential of itinerant 
particles increases. 
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APPENDIX: EXACT RELATIONS BETWEEN 
GREEN'S FUNCTIONS 

Let us consider the effective single-site problem in 
terms of the auxiliary Fermi-field 

m 

H cS =H + H i + J2v a (aU„ + ( A. 1 ) 

(7=1 

Hq is a single-site Hamiltonian; is an auxiliary envi- 
ronment Hamiltonian. The number m of sorts of itiner- 
ant particles can be arbitrary (m = 1, 2, 3, . . .). So, the 
effective Hamiltonian can describe the Falicov-Kimball 
model (m = 1 ) and the Hubbard model (m = 2). The 
algebra of £ operators is defined by the anticommutation 
relations 

= = °> M-y=^. (A.2) 

It can be proved that the following relations take place 

V a {(Z„\A)) u = Mu)({a„\A)) u , (A.3) 
27rt£«k|£)) u = Ja{ul ) + 2 7 rj2H«a>t)) L ,,(A.4) 

where A is the arbitrary Fermi- operator that anticom- 
mutates with the £ operators, and J CT (w) is the Green's 
function for the unperturbed Hamiltonian 

Proof. Thermodynamic perturbation theory can be 
formulated based on the interaction representation for 
the statistical operator with the use of the Hq + H^ op- 
erator as the zero-order Hamiltonian: 



P : 



-PH oil _ -P(H +Ht:) 



HP), 



dTH int (r) 



(A.5) 
(A.6) 



The part of the interaction Hamiltonian describing one 
sort (tr) of particles is separated off: 



(A.7) 
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The residue B a commutates with the operators of the 
chosen sort a: [B a , £ CT ] — [B a , = 0. Thus, at the 

perturbation expansion the operator does not have 
to be paired with the operator B a . Here we introduce 
notations for the Green's functions 

Mr-r') =(T t &(t)Ut j ))o 
G a (r-r')=(T T aUr)a a (r')) 

Q„(t-t j ) ={T T A a {T)Wr')) 
G*(T-T>) = (T T AjT)ajT<)) | 



The perturbation theory expansion of the scattering ma- 
trix a((3) gives the following series for the Green's func- 
tion $ CT (r - t')\ 



V? +1 ^ (-1) ^ (2p+l)! f p , , 
' W/3))o^(2 P + l)!^/!(/ + l)!(2p-20!7o 



x (T T A(T)Uy)(al^), ■ ■ ■ ■ • ■ {£a a ) 2l+l B a {T 2l+2 ) ■ ■ ■ B a (r 2p+1 )) . (A.8) 

The averaging of T products is performed in the zero-order Hamiltonian according to the Wick's theorem by the 
consecutive pairing. We start the pairing procedure from the operator ^(r') and it is performed only with the 
operators £j. After the first pairing we have the following expression: 

y2l + l r!3 aa p ^ „p 

^ {T - T ' ) = vw,L dr ^ (T ''" T \?oS^^^ dri -- dT2p 

x {T T A{r)a rJ {r"){alQ 1 ■ ■ ■ (aU a )M a *)l+i ' ' ' (#0 sA^+O ' ' ' S CT (r 2p )) . (A.9) 

I 



Summing up the series in i|A.9|) gives the Green's function 
G^(t — t') and we have 

$ ff ( T - r ') = F CT f 1 dr"Gi(T-T")g a (T" '-r'), (A.10) 
Jo 

or in the Matsubara frequency representation: 

= VaGp (id n )Qo-(u} n ) . (A.ll) 



Finally, in order to obtain the expression l|A.3() . an 
analytical continuation of the Green's functions from 
the imaginary axis to the real one should be done 
(ioj n -> w + ie, G^{oj n ) -> 27r((a CT |l)) w , $ CT (w„) -> 
27rV r CT ((^|l)) w , ...)• If we put A = Q, the pairing of 
the operator ^ CT (r') with £J.( r ) gives the first term in 
(|A.4(I . The second term is obtained using the procedure 
the same as at deriving the relation l|A.3|) . 
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